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We consider the linear stability of two-dimensional nonlinear magnetohydrody- 
namic basic states to long-wavelength three-dimensional perturbations. Following 
Hughes & Proctor (2009a), the 2D basic states are obtained from a specific forcing 
function in the presence of an initially uniform mean field of strength B. By ex- 
tending to the nonlinear regime the kinematic analysis of Roberts (1970), we show 
that it is possible to predict the growth rate of these perturbations by applying 
mean field theory to both the momentum and the induction equations. If S = 0, 
these equations decouple and large-scale magnetic and velocity perturbations may 
grow via the kinematic a-effect and the AKA instability respectively. However, if 
B 0, the momentum and induction equations are coupled by the Lorentz force; in 
this case, we show that four transport tensors are now necessary to determine the 
growth rate of the perturbations. We illustrate these situations by numerical exam- 
ples; in particular, we show that a mean field description of the nonlinear regime 
based solely on a quenched a coefficient is incorrect. 
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1. Introduction 

Astronomical observations show that most astrophysical objects are magnetised and 
often possess a coherent large-scale magnetic field. Such fields are believed to be 
the results of a large-scale dynamo, whereby inductive motions within a conductive 
fluid can generate and sustain a magnetic field on scales much larger than their 
own. 

The first step in addressing this problem is to understand how a weak seed field 
can be amplified by the motion. To study this kinematic phase, it is assumed that 
the magnetic field has no effect on the motion, which can therefore be considered 
as prescribed. Kinematic large-scale dynamos are traditionally studied within the 
framework of mean field electrodynamics, a turbulence closure theory in which the 
effects of small scale interactions are parametrised by transport coefficients. Within 
this theory, the equation for the evolution of a large-scale magnetic field B driven 
by a motion U (with no mean) is typically of the form 

dtB = V X {a-B + 13 -VB + Rm-^V xB) , (1.1) 

where a and (3 are transport pseudo-tensors that depend on the characteristics of 
U and on the magnetic Reynolds number Rm (sec, for example, Moffatt 1978). The 
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symmetric part of a is the so-called a-effect, which is responsible for the growth 
of B. For isotropic turbulence, we have simply that = aSij and (iijk = I3^ijk\ 
in this case, the growth rate cr of a long-wavelength mean field with wavenumber 
A; <C 1 can be approximated by 

a Pi ak - {P + Rm-^)k^ . (1.2) 

(For completeness, it should be noted, as discussed by Hughes & Proctor (2009b), 
that although equations (1.1) and (1.2) are entirely consistent, there is a potential 
inconsistency between the diffusion term in (1.2) and that derived from considering 
the decay of a large-scale field. This, however, is not germane to the problem studied 
in this paper, in which we consider only effects of first order in k.) 

If cr > 0, an initially weak large-scale magnetic field will grow exponentially 
until a stage is reached at which the back-reaction of the magnetic field on the flow 
U must be taken into account. This implies that the induction equation has to be 
solved in conjunction with the momentum equation, in which the Lorentz forces 
are included. The hope, however, has been that mean field electrodynamics is still 
applicable with B still satisfying (1.1), but with transport coefficients that now 
depend on the mean field, i.e. ot = ot{B) and /3 = (3{B). Over the past two decades, 
considerable effort has been devoted to the analytical and numerical determination 
of these coefficients and their exact dependence on B, the so-called a-quenching 
issue (e.g. Vainshtein & Cattaneo 1992, Kulsrud & Anderson 1992, Gruzinov & 
Diamond 1994, Cattaneo & Hughes 1996, Blackman & Field 2000; an extensive list 
of references may be found in the review by Brandenburg & Subramanian 2005). 
Recently, it has been suggested that the dependence of o; = cn{B) and /3 = l3{B) 
on B can be determined by considering only the effects of the velocity field modified 
by the Lorentz force (Tilgner & Brandenburg 2008). In this paper we argue that 
such an approach, and indeed any approach that privileges the magnetic field over 
the velocity field in a magnetohydrodynamic (MHD) state, is flawed. We explain 
how mean field electrodynamics can be applied to the nonlinear regime only if the 
coupling between the induction and the momentum equation via the Lorentz force 
is fully taken into account. As a result, new transport tensors emerge that are 
simply not captured by considering the induction equation alone. 

Our aim therefore is to consider long-wavelength perturbations to a nonlinear 
MHD state, and to explain how the growth (or decay) of such perturbations can 
be explained in terms of transport coefficients arising from small-scale interactions. 
This can be regarded as the dynamic analogue of the kinematic dynamo problem, 
in which the growth of a long-wavelength magnetic perturbation to a purely hy- 
drodynamic state can be described in terms of the a-tensor. The crucial feature of 
our analysis of a coupled MHD state is that the velocity and magnetic fields must 
be treated on an equal footing throughout; failure to do this will lead to physically 
meaningless results. 

In order to make clear the essential features involved, we shall in this paper 
confine our attention to the restricted problem of investigating 3D long- wavelength 
linear perturbations to 2D basic MHD states — the basic state variables depend on 
X, y and t, the perturbations are taken to have a long wavelength in the z direction. 
There are three great advantages in considering this slightly restricted 2D problem. 
The first is that it is conceptually simpler than considering fully 3D basic states; 
the second is that average quantities are readily calculated; the third is that it 
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allows a direct comparison between the growth rates determined via evaluation 
of the various transport tensors and those calculated by a direct solution of the 
stability problem. Although the generalisation to 3D is formally straightforward 
(the corresponding equations are derived in Appendix A), there are a number of 
subtle issues, not present in 2D, that arise in the measurement and interpretation of 
the transport tensors. We shall discuss these in a future paper (Courvoisier, Hughes 
& Proctor 2010). 

We follow the approach recently taken by Hughes & Proctor (2009a), who con- 
sidered the growth of linear 3D perturbations of 2D MHD basic states obtained 
from a specified forcing function in the presence of an imposed initially uniform 
magnetic field of strength B. Crucially, their analysis involved solving a system of 
coupled equations for both magnetic and velocity perturbations. In this paper, we 
use exactly the same framework but, this time, we are interested in the stability 
properties of long- wavelength perturbations. Although these are not in general the 
dominant modes of instability, they arc the modes that can be rigorously described 
using mean field theory, as discussed above. If the imposed field is small (S <C 1), 
the velocity field of the basic state, U say, is determined solely by the prescribed 
forcing. We then recover the kinematic dynamo problem, for which mean field elec- 
trodynamics can be readily applied and in which large-scale magnetic perturbations 
can grow via the a-effect, provided U lacks parity invariance. Furthermore, if U 
is also anisotropic, large-scale velocity perturbations may grow via the anisotropic 
kinetic alpha (AKA) instability (Frisch, She & Sulem 1987). For larger values of 
B however, the induction and momentum equations are coupled and cannot be 
considered separately. For this case, we develop a 'nonlinear' theory of mean field 
electrodynamics . 

The format of the paper is as follows. In § 2, in order to investigate, in gen- 
eral terms, the long-wavelength 3D linear instabilities of 2D nonlinear MHD basic 
states, we pursue an analysis d la Roberts (1970), who considered the a-effect in 
the kinematic regime. The main result of our paper — the derivation of the four 
transport tensors necessary for the description of long- wavelength perturbations — 
is contained in § 2(b). In § 2(c) we show how, for a purely hydrodynamic basic state 
(no imposed magnetic field), the evolutions of the magnetic and velocity pertur- 
bations decouple, leading to the possibility of a kinematic dynamo (for magnetic 
perturbations) or the AKA instability (for velocity perturbations). In § 2(d) we dis- 
cuss what might be thought of as the kinematic dynamo problem for a fully coupled 
MHD flow — in other words, considering magnetic perturbations but, arbitrarily, 
neglecting velocity perturbations. Although this problem has received some atten- 
tion, we argue that such an approach is physically inconsistent. We then present 
two numerical examples. The first, described in § 3, illustrates our theory for both 
hydrodynamic and magnetohydrodynamic basic states. The second, detailed in § 4, 
is designed to emphasise that considering the a-effect alone fails to describe accu- 
rately large-scale dynamo action in the nonlinear regime. Finally, we discuss the 
implications of our results in § 5. 
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2. Analytical approach 

(a) The 2D basic state 

We consider, with the usual notation, the basic state B{x,y,t), U{x,y,t), 
P{x, y, t) resulting from the solution of the incompressible MHD equations, 

{dt- Rm-^V"^) B = V X {U X B) + B-VU, (2.1) 
{dt - i?e-i V^) U = -VP + V • {BB - UU) + B VB + F, (2.2) 
V-B = 0, V-C/ = 0, (2.3) 

where the forcing F{x,y,t) is spatially periodic and B is a constant vector field; 
Re and Rm are, respectively, the kinetic and magnetic Reynolds numbers. Since 
dynamo action cannot regenerate a 2D magnetic field (Cowling 1933), B is due 
solely to the distortion of B by the fluid motion. Furthermore, since U is indepen- 
dent of z, the mean flux B needs to be horizontal in order to have a non-zero B 
(for convenience, we shall refer to the xy-plane as the horizontal plane, although 
gravity plays no part in the present study). 

In the analysis that follows, we consider basic states that are both spatially and 
temporally periodic, for which it is natural to define an averaging operation by 

(U) = J^l I I U{x,y,t)dxdydt, (2.4) 

where the flow has spatial periodicity 2tt and temporal period T. We restrict our 
analysis to the case where (B) = (U) = 0, i.e. there is no mean flow and the mean 
magnetic field corresponds to B. 

Of course, systems described by (2.1) - (2.3) can exhibit complex temporal 
behaviour even with a time-periodic forcing, as shown by Hughes & Proctor (2009a). 
In such cases, our analysis is still valid provided that averages arc taken over times 
long compared with the characteristic time scale of the fiuctuations in the basic 
state. 

(6) Growth rate of z- dependent perturbations 

We now consider small-scale, incompressible, linear perturbations to the basic 
state, of the form 

{b{x,t),u{x,t),p{x,t)) = (if(,x,y,t), V(,x,?y,t),n(.T,y,t)) e^'=^+''W*, (2.5) 

where V, H and 11 vary on the same spatial and temporal scales as the basic state, 
and where p{k) = a{k) + iuj{k) is the complex growth rate of the perturbation, 
which depends on its wavcnumbcr k. 

The perturbations can be decomposed into average and fiuctuating parts; fol- 
lowing Roberts (1970) we define the average by 

(b{z,t)M^.,t),p{z,t)) = {{H),{V),{U)) e^'^^+fW*. (2.6) 

Our aim is to illustrate how a mean field approach can be used to determine the 
growth rate p{k) of modes with non-zero mean magnetic fields, i.e. modes for which 
{H) ^ 0, in the hmit of A; ^ 0. 
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The equations satisfied by the linear perturbations are given by 

(p + Rm-^k^) H+{dt- Rm-^V^) H = {W + ikz ) x {V x B + U x H) 

+ B-WV, (2.7) 
{p + Re-^k^) V+{dt- Re-^V^) F = (V + zA:i ) • {HB + BH -VU - UV) 

-Vn-ifcHz +B VH, (2.8) 

V ■ H + ikz ■ H = Q, V -V + ikz - V = Q. (2.9) 

We now treat fc as a small parameter and expand H, V, 11 and the growth rate p 
in powers of fc, so that, for example, 

H = Ho + H^ + ...Hn + ..., (2.10) 

where H „ is of n*^ order in k. Each H,, can be decomposed into its mean and 
fluctuating parts, i?„ = {Hn) + h^, similarly, Vn = (^n) + Vn- 
At zeroth order, equations (2.7) - (2.9) become 

po-ffo + {dt - Rm-'^V^) Ho^W X {Vox B + U X Ho) + B- VVq, (2.11) 
PoVo + {dt - Re-^V^) Fo = V • {H^B + BH^ - V^U - UVq) 

-Vno + B-V/To, (2.12) 
V-Ho = Q, V-Vo = 0. (2.13) 

On averaging these equations, we obtain 

po(-H'o) =Po(Vo) =0. (2.14) 

For non-zero mean field solutions it follows therefore that po = 0. The corresponding 
equations for the fluctuations then become 

(dt - Rm-'^V^)ho = Vx{voxB + Uxho) 

+ {Ho)-VU -{Vo)-VB + B-S7vo, (2.15) 
{dt - -Re" V^)vo = -VHo + V • {hoB + Bho - vqU - Uvq) 

+ (ifo) • VB - (Vo) • VJ7 + B • V/io, (2.16) 

V ■ho = 0, V • -uo = 0. (2.17) 

The first non-trivial mean field equations are obtained from the horizontal averages 
of equations (2.7) - (2.9) at order k; we obtain 

pi{Ho) =ikz X {vox B + U xho), (2.18) 
pi{Vo) = -iHokz +ikz ■ {hoB + Bho - vqU - Uvo), (2.19) 
z ■ {Ho) =0, z ■ {Vo) = 0. (2.20) 

It can be seen from equations (2.15) - (2.17) that ho and Vq are subject to a Hnear 
forcing by both (Hq) and {Vq). Therefore, (2.18) and (2.19) can be rewritten 
(without approximation) as 

Pi{Ho) = ikz X (a^ • {Ho) + • {Vo)) , (2.21) 

Pi{Vo) = ~ikz Ho + ikz ■ (r^ • {Ho) + ■ (Fo)) , (2.22) 
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or, in component form, as 

Pi{Ho)i = -ikeaem «„(^fo)n + a^JVo)n) , (2.23) 
Pi{Vo)i = -ikSesIlo + ik (Tie^{Ho)m + T^im{Vo)m) . (2.24) 

The tensorsf cx^ , , and are constant and depend on the basic state U 
and B — hence on the forcing F, the imposed magnetic field B and the Reynolds 
numbers. The tensors and are symmetric with respect to their first two 
indices. 

In order to use mean field theory to determine the growth rate of long- wavelength 
perturbations to an MHD state, it is necessary to determine the four tensors oF , 

, and by solving equations (2.15) - (2.17) for Hq and Vq, with imposed 
constant magnetic and velocity fields, alongside the basic state equations (2.1) - 
(2.3). 

Wc not(^ that for the 2D system considered here, the generated mean fields are 
horizontal. Furthermore, only the horizontal component of the mean emf, {vq x 
B + U X ho), contributes to the growth of the mean magnetic field; therefore only 
the 2x2 parts of a'^ and a.^ relating horizontal quantities need be calculated. 
Similarly, only the horizontal component of the right hand side of (2.19) influences 
pi, so the only components of the mean stress tensor of interest are of the form 
{ho3 Bf + Bs hoe — '-'03 Ue — C/3 voe), with I = 1,2. Thus we need only to determine 
^Mm = ^I3m and Vf,^ = Vg^ for i,m=l, 2. 

Equations (2.21) and (2.22) can then be rewritten as 
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(2.25) 



The first order growth rate is then an eigenvalue of the matrix Ad k. 

In §§ 3 and 4, we test our analysis numerically by comparing the actual growth 
rate of .z-dependent perturbations with predictions obtained by calculating the co- 
efficients of the four tensors oc^ , oc^ , and . Prior to this though, we consider 
two special cases of the above theory. 



(c) Purely hydrodynarnic basic states, 3=0 

In this subsection, we consider the limit when B = 0. This implies that B = 
and the basic state is purely hydrodynamic. In this case, equations (2.15) and (2.16) 
for the magnetic and velocity perturbations simply become 

{dt - Rm-^V^)ho =\7 x{U xha) + (Ha) ■ VU, (2.26) 
[dt - Re-^V^)vo = -VHo - V • {voU + Uvo) - (Fo) • VC/. (2.27) 

t To be more precise, c»F and are true tensors, whereas and are pseudo-tensors. 
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The corresponding first order equations giving the growth rate, (2.18) and (2.19), 
become 

pf{Ho)=ikz x{Uxha)=ikz x{a.-{Ha)), (2.28) 
p\{Vq) = -iUok z -ikz ■ {voU + Uvo) = ikz ■ {T ■ {Vq)) . (2.29) 

Wc sec that in this Umit the evohitions of ffo and Vo decouple, and so wc expect 
the growth rates pf and to be distinct. The magnetic problem, described by 
(2.26) and (2.28), reduces to the kinematic dynamo problem for the velocity field 
U , with OL the standard a. tensor of mean-field electrodynamics. The symmetric 
part of this (pseudo)-tensor encompasses the a-effect of kinematic dynamo theory; 
it can be non-zero only in flows lacking parity invariance, such as helical flows (see 
Moffatt 1978). The hydrodynamic problem, described by (2.27) and (2.29), reduces 
to that studied by Frisch et al. (1987). The coefficient F here corresponds to the 
anisotropic kinetic alpha (AKA) effect, which can be nonzero in anisotropic flows 
lacking parity invariance. This effect also requires that a uniform large-scale velocity 
field produces a change in the average small-scale Reynolds stresses, therefore that 
Galilean invariance be broken. For more details and a comparison between the 
conditions required for the existence of the a and AKA effects, see Frisch et al. 
(1987). 

Equations (2.28) and (2.29) can be rewritten in matrix form as 









p? 


Hoy _ 






' Vox ' 
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(2.30) 

(2.31) 



so the growth rates pf and pY are eigenvalues of the matrices A.k and Qk, re- 
spectively. The a-cffcct can be readily determined by solving (2.26) for an imposed 
horizontal mean field (Ho), while T can be obtained by solving (2.27) for an im- 
posed horizontal mean flow (Vq). 

The limiting procedure under consideration here (i.e. taking 3 = 0) constitutes 
a well-defined physical problem. It corresponds to the case when there is no Lorentz 
force (i.e. in 2D there is no mean flux, or, in 3D, Rm and Re are low enough 
to prevent small-scale dynamo action and small-scale hydrodynamic instabilities); 
the momentum and induction equations therefore decouple and the a and AKA 
instabilities can exist independently. A numerical example of a forcing for which 
both a and AKA-cffects are present is given in § 3, in which we also describe what 
happens when 0^0. 



(rf) Magnetic perturbations only, u = 

In this subsection, we consider the problem of perturbing the basic state only in 
the magnetic field, i.e. we arbitrarily set u = 0, and hence = 0. This corresponds 
to studying a kinematic dynamo problem for the basic state velocity field U{B). 
This line of research has been considered recently by Cattaneo & Tobias (2009) and 
by Tilgner & Brandenburg (2008). 
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In this case, equation (2.15) for the magnetic perturbation at zeroth order simply 
becomes 

{dt - Rm-^W'^) ho = V X {U{B) x ho) + (Hq) ■ VU{B). (2.32) 

Our notation here, which shows that the basic state flow U, and hence ho, depends 

on S, is chosen to emphases the difference between (2.26) and (2.32) 
Equations (2.18) and (2.21), which give the growth rate, become 

(iJo) = ikz X {U{B) X ho) ikz x {a{B) ■ (Hq)) , (2.33) 

where a{B) corresponds to the a tensor of mean field electrodynamics. It depends 
on the mean flux of the basic state, B, as well as on F and the Reynolds numbers; 
in that sense, it is distinct from the purely kinematic a tensor described in §2(c). 
We expect ot{B) to be different from cx^ in general, since the latter encompasses 
the dependence of both ho and vo on the mean magnetic field {Ho). 

Using (2.33) we can define the matrix .4.(6) whose eigenvalues give the growth 
rate pf. Its entries are similar to those of A., but again we emphases the B depen- 
dence in the notation. 

The cx{B) tensor can be readily determined by solving equation (2.32). We note 
that this equation is formally equivalent to the induction equation (2.1) from the 
basic state, with ho = B and {Hq) = S; therefore some components of a{B) can 
be directly computed from the mean emf {U{B) x B) emerging from the basic 
state. To obtain the other components of a, it is necessary to solve a 'passive' 
induction equation alongside the basic state, with the imposed mean field in a 
direction perpendicular to that of S, i.e. in the present context, this corresponds to 
solving equation (2.32) with (Hq) = {Ho) y ■ This procedure for calculating a{B) 
is now often referred to as the 'test-field' method (Schrinner et al. 2005) and has 
been discussed in the context of a-quenching by, for example, Brandenburg et al. 
(2008), Tilgner & Brandenburg (2008). 

It is important to note that we can use the basic state to determine the kinematic 
a-effect for U{B) only because the system is two-dimensional. Indeed, equation (2.1) 
is part of a system of coupled equations, while (2.32) is a passive vector equation 
for the same velocity field. The 3D simulations of Cattaneo & Tobias (2009) show 
that, despite being identical, these equations in general have different magnetic 
field solutions. In 2D however, we find that both equations converge to the same 
solution, owing to the fact that dynamo action is not possible in this case and that 
therefore the solution is driven solely by the imposed mean magnetic field. 

Finally, we would like to stress the important point that considering only mag- 
netic perturbations constitutes a problem with no physical meaning, since, as dis- 
cussed in the introduction, it privileges the magnetic field over the velocity field in 
a fully coupled MHD state. Crucially, there is no physical justification for setting 
V = 0. Indeed, from (2.16), this would imply that V • {hoB + Bho) + B ■ Vho = 0, 
which, in general, is in contradiction with the assumed existence of non-zero mag- 
netic perturbations. 

As argued by Hughes & Proctor (2009a), it is therefore necessary to consider 
perturbations to both the induction and momentum equations. In this case, the 
actual growth rate of the fiuctuations is completely unrelated to the kinematic a- 
eSect corresponding to U{B), as we shall demonstrate numerically in §4. Instead, 
all four tensors , a^, and need to be determined. 
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B 


relevant matrix 


eigenvalues 





A 

Q 


±0.077 

±0.386 


0.1 


M 


±0.209 + 0.182 i ±0.171- 0.176 i 



Table 1. Eigenvalues of the matrices for the AKA forcing; Re = 1, Rm = 16. 

For c;oniplctcness it should perhaps bo noted that there is a related iinphysical 
problem, which would consist of investigating the stability of velocity perturbations 
to an MHD state, disregarding magnetic perturbations. This would involve the 
solution of (2.16) but with ho = Hq = 0, i.e. 

{dt - Re-^V^)vo = -VHo - V • {voU{B) + U{B)vo) - {Vo) ■ VU{B). (2.34) 

3. Numerical example I: AKA forcing 

(a) Basic state and mean field coefficients 
For this example, we choose the force in the momentum equation (2.2) to be 

F= (jie-'^ \/2 cos {y + Re-H),Re-'^ \/2 cos {x - Re-^t),F^ + (3.1) 

This forcing was first proposed by Frisch et al. (1987) in order to drive a flow 
displaying the AKA-cfFect. By design, flows driven by this forcing lack parity in- 
variance, and thus, in general, may be expected to drive both a and AKA effects. 
Frisch et al. (1987) determined the F tensor analytically when i?e <C 1, thus allow- 
ing the neglect of the inertia terms from the momentum equation. Here we drive 
the flow with Re = 1; in this case, and in contrast to when i?e <C 1, the helicity 
of the ensuing flow is non-zero. We flrst consider the case when B = 0, described 
in §2(c), and restrict our simulations to Rm, = 16. In this parameter regime, both 
U and B are time-periodic with a well-defined frequency. The matrices A and Q 
are determined by solving (2.26) and (2.27) respectively for solenoidal vector fields, 
which is equivalent to solving the full system (2.15) (2.17), since B = 0. The 
matrices are given in Appendix B and their eigenvalues are contained in table 1. 
These are real and show that the velocity field considered does indeed display both 
a and AKA effects; the latter is of greater magnitude, which is maybe not unex- 
pected for a velocity field specifically designed to drive an AKA effect and with 
little kinetic helicity, which, although, not strictly necessary for an a-effect, is often 
helpful (Gilbert, Frisch & Pouquet 1988). 

We then go on to consider the case when B = 0.1x , which leads to the coupling 
of magnetic and velocity perturbations via the Lorentz force. We determine the 
matrix Al, the elements of which are given in Appendix B, and its eigenvalues, 
which are given in table 1; they show that the growing modes are now oscillatory. 

(6) Growth of z- dependent perturbations 

For comparison with the mean field theory results, we calculate directly the 
growth rate of z-dependent linear perturbations to the basic state. Following the 
method of Hughes & Proctor (2009a), we solve equations (2.7) - (2.9) numeri- 
cally, using a time stepping procedure, in order to determine H{x, y, t) e^^'^'* and 
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Figure 1. Results for the AKA forcing; B = Q, Re = 1, Rm = 16. Left: Time series of the 
energy in the perturbed magnetic and velocity fields for k = 0.004. The dashed lines show 
the predicted growth rate for the a and AKA instabilities. Right: crik) for the magnetic 
(+) and velocity (★) perturbations. The dashed lines correspond to A k, where A stands 
for the positive eigenvalue of A. (for the magnetic perturbations) or of Q (for the velocity 
perturbations). 



V{x, y, t) e^'('^)*. We obtain a{k) = Re\p{k)] from a least squares fit to the time series 
for the corresponding magnetic and kinetic energies, and determine the frequency 
u){k) = Im[p{k)] by inspection of the time evolution of spacc-avcragcd magnetic; or 
velocity components. Since we are interested in verifying our predictions for p{k) 
at low k, we only consider values of k less than 0.1. 

The; rc;sults for 6 = arc given in figure 1. As discussed above, the induction and 
momentum equations decouple, thus allowing magnetic and velocity perturbations 
to evolve independently. The left panel of figure 1, which presents typical time series 
of the energy in the perturbed magnetic and velocity fields (here for k — 0.004), 
confirms that these modes have different growth rates. The right panel summarises 
the results obtained for a{k) and shows that, for low values of k, a Xk, where A 
stands for the positive eigenvalue oi A.or Q, for magnetic or velocity perturbations 
respectively, in agreement with the theory derived in §2(c). 

The results for B = 0.1 are given in figure 2, which presents a{k) (left panel) 
and uj{k) (right panel) together with _Re[A] k and /m[A] k (dashed lines), where A 
is the eigenvalue of Ad with the largest real part. Here again we see that the linear 
behaviour of both a{k) and ui{k) at low k is well predicted by the analysis. 

Finally, wc note that the dynamo problem, including the nonlinear saturation of 
the instability, in a flow driven by the AKA forcing (or similar) has been considered 
by, for example, Galanti, Sulem & Gilbert (1991) and Brandenburg & von Rekowski 
(2001). In our notation, this would correspond to studying the saturation of the 
instability in the case when B = 0, i.e. when the basic state is purely hydrodynamic. 
The focus of the present paper is however the linear phase of the instability when 
B^O. 
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Rm=16, Bo = 0.1 Rm=16, B„ = 0.1 




Figure 2. a{k) (left) and w(fc) (right) for the AKA forcing; Rm = 16, B = 0.1. The dashed 
lines corresponds to Re[X] k (left) and Im[\]k (right), where A is the eigenvalue of M. 
with the largest real part. 

4. Numerical example II: MW+ forcing 

(a) The basic state 

We now choose F so that, in the absence of B, it drives the MW+ flow of Otani 
(1993), given by 

Uo{x,y,t) = {dyi:,-d^i^,-^P), (4.1) 

with 

ip{x, y, t) = cos X cos^ t — cos y sin^ t. (4.2) 

This flow has the Beltrami property that U is parallel to V x t/ and so F is simply 
the solution of 

F = dtUo - Re-'^V'^Uo. (4.3) 

To ensure that the flow driven when B = settles down rapidly to Uq and to 
avoid any hydrodynamic instabiHties we restrict our computations to Re = 1. The 
symmetry properties of the MW+ flow imply that there is no AKA effect and that 
its kinematic a tensor is diagonal, with entries an = (X22 = a; these have been 
computed numerically by Courvoisier (2008). 

We use this example to illustrate the fact that taking only magnetic perturba- 
tions into account can be misleading. To this end, we consider both the full problem 
and the artificial case of taking V = 0, discussed in § 2(d). As in the previous sec- 
tion, we take B in the a;-direction. For the parameters considered here, the basic 
state proves to be time-periodic with a well-defined frequency. 

We first compute the eigenvalues of Jk{B) by using the basic state and solv- 
ing (2.32) with an imposed mean field in the direction perpendicular to B. The 
results are given in table 2 (see Appendix B for the full matrices). For B <C 1, 
A{B) ^ A., and we recover the kinematic a-effect for the MW+ flow. As B in- 
creases for Rm = 16, we see that the eigenvalues of A.{B) decrease in magnitude, 
as expected from a-quenching studies. This tells us that the flow U{B), taken in 
isolation, is not as good a kinematic large-scale dynamo as the flow U. 
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Rm 


B 


eigenvalues 


16 


0.001 


±0.988 




0.25 


±0.624 




1.0 


±0.088 


128 


0.1 


±1.110 



Table 2. Eigenvalues of A{B) for the MW+ forcing; Re = 1. 



Rm 


B 


eigenvalues 


eigenvectors 


16 


0.001 


±0.988 



(Ho) only 




0.25 


±0.085 
±0.406 


(Ho) only 
(Vo) only 




1.0 


±0.217 i 

±0.437 


(Ho) only 

{Vo} only 


128 


0.1 


±n.:M2 

= 1.27() / 


(Ho) only 
(y,j) only 



Table 3. Eigenvalues of Ai for the MW+ forcing; Re= 1. 

We then consider the full problem, as described in §2(b), and compute the 
eigenvalues of by solving the system (2.15) - (2.17); the elements of Ai are 
contained in Appendix B. It appears that for all the values of B and Rm investigated 
here, its eigenvectors arc cither purely magnetic or purely hydrodynamic, so that 
the mean magnetic and velocity fields do not grow at the same rate. The results 
for the eigenvalues and the nature of the corresponding eigenvectors are given in 
table 3. For Rm = 16, B = 0.25 and 1, the eigenvalue with the largest real part 
corresponds to a mean flow solution. This might be due to the fact that a non-zero 
mean flux introduces the anisotropy favourable to the development of an AKA-type 
instability. For B = 0.25, the eigenvalue corresponding to a mean magnetic field 
solution is pure imaginary, so this mode is marginally stable here, in contrast to 
the growing mode found for magnetic only solutions. For Rm = 128 and B = 0.1, 
the situation is reversed, as the eigenvalue with the largest real part corresponds 
to a mean magnetic field solution, possibly because higher values of Rm (for fixed 
Re) favour magnetic field growth. 

(6) Growth of z- dependent perturbations 

We now calculate the growth rate of z-dcpendcnt perturbations. Wc consider the 
artificial magnetic problem, for which we impose V = 0, as well as the full problem, 
in which magnetic and velocity perturbation equations are solved in concert, as in 
Hughes & Proctor (2009a). For the latter wc arc able to determine the growth rate 
of both the magnetic and velocity eigenvectors by following the time evolution of 
the spatially averaged H{x, y, t) e^^'^^* and V{x, y, t) e^'^'^^*. 

Figure 3 shows cr{k) for Rm ~ 16, with B = 0.001 (left panel) and 0.25 (right 
panel). The + correspond to the fastest growing mode in the magnetic problem, 
while the ★ and o correspond to the dominant velocity and magnetic modes, re- 
spectively, in the full problem. Figure 4 shows the results for Rm = 16 with B = 1 
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Rm=16, B„ = 10" 



0.10 



0.04 



0.02 



0.00 



0.00 0.02 0.04 0.06 0.08 0.10 
k 



Rm=16, B„ = 0.25 




0.00 0.02 0.04 0.06 0.08 0.10 
k 



Figure 3. Left: cr(fc) for Rm = 16, B = 0.001. The dashed line corresponds to \ k, where 
A is the positive eigenvalue of A.. Right: cr(fc) for the magnetic problem (+) and the full 

problem (★, velocity mode, and o, magnetic mode) for Rm = 16 and B — 0.25. The dashed 
lines corresponds to A k, where A stands for the positive eigenvalue(s) of A.{B) (for the 
magnetic problem) and of M. (for the full problem) . 



and for Rm = 128 with B = 0.1. The left panels show a{k) for the fastest growing 
mode in the magnetic (+) and full (★) problems. The right panels show u}{k) for 
the marginal mode in the full problem. In all cases, the dashed lines give the rele- 
vant linear approximation for low k based on the eigenvalues of Jl{B) or Ai. This 
shows that our results are in excellent agreement with the analysis presented in 
§ 2. These results also demonstrate that in the cases investigated here, the growth 
rates of magnetic perturbations taken on their own are very different from those 
of the combined magnetic and velocity perturbations, which we believe to be the 
relevant situation to study. As shown by this particular example, this difference 
persists even in situations for which the relevant components of a'^ and vanish 
(see Appendix B). 



5. Conclusions 

In this paper wc have investigated the instability of two-dimensional magneto- 
hydrodynamic states to long-wavelength three-dimensional disturbances. As well 
as calculating the growth rates directly, we have used the methods pioneered by 
Roberts (1970) for the simpler kinematic dynamo, in order to derive a mean field 
theory for the linear growth of these perturbations, with coefficients that can be 
determined by imposing spatially uniform perturbations on the basic state. The 
calculated transport coefficients in the mean field theory give excellent agreement 
with the directly computed results. 

The crucial feature of the mean field equations is that perturbations in the veloc- 
ity and magnetic field are coupled, and that, in general, the growing cigcnfunctions 
involve both sorts of field. Thus it is imperative that the induction and momen- 
tum equations are considered on the same footing. In contrast, earlier studies of 
the effects of equilibrated fields on the traditional a-effect have used the induction 
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Rm=16, B„ = 1 Rm=16, B„ = 1 




0.00 0.02 0.04 0.06 0.08 0.10 0.00 0.02 0.04 0.06 0.08 0.10 

k k 



Figure 4. Results for the MW+ forcing; Rm = 16, B = 1.0 (top plots); Rm = 128, B = 0.1 
(bottom plots). Left: (j{k) for the magnetic problem (+) and the full problem (★). The 
dashed lines corresponds to f?e[A] k, where A stands for the eigenvalue with the largest real 
part of A.{B) (for the magnetic problem) and of Al (for the full problem). Right: u!{k) 
for the full problem. The dashed lines corresponds to 7m [A] k, where A is the eigenvalue 
of M. corresponding to the marginal mode. 

equation without any coupling to the momentum equation, and this 'passive' theory 
gives incorrect results in general. 

Our results shed a new Ught on the general problem of the long-wavclcngth 
instabilities of saturated MHD states. In the past, the sole focus of interest has 
been the nature of the a-effect, and its dependence on ambient magnetic field 
levels (a-quenching). However, we can now sec that there are situations in which 
the preferred unstable eigenmode may be dominated by the velocity field, and, 
furthermore, that there are cases for which magnetic and velocity perturbations 
grow faster than magnetic perturbations taken on their own. Thus even when strong 
('catastrophic') quenching is predicted on the basis of the old ideas, we now see that 
there could be other mechanisms leading to a large-scale instability. 

Here we have concentrated on 2D examples, which, although somewhat artificial. 
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have the great advantage of aUowing us to compare the predictions of mean field 
theory with the actual growth rate of the magnetic field, because we are able to 
isolate a single mode artificially something that would be much harder to achieve 
in 3D with the current computational resources. However, we emphasize that our 
main findings apply to the linear long-wavelength instability of any nonlinearly 
saturated MHD system, whether in two or three dimensions. There are however 
important differences between the 2D and 3D derivation of the formal 

mean field equations in the 3D case is given in Appendix A, and a full treatment 
will be given in a subsequent paper (Courvoisier et al. 2010). 



We are extremely grateful to Andrew Gilbert and Steve Tobias for many enlighten- 
ing discussions. This work was supported by STFC and by a Royal Society Lever- 
hulme Trust Senior Research Fellowship (DWH). Some of the work was performed 
during the 2009 workshop on Dynamo Theory at the Institut Henri Poincare; we 
are very grateful to the organisers, Emmanuel Dormy, Stephan Fauve and Frangois 
Petrelis. 



Appendix A. The general case of a 3D basic state 

Here we derive the more general equations for a 3D basic state B{x,t), U{x,t), 
P{x,t), which we shall assume, for simplicity, to be spatially and temporally pe- 
riodic. The calculations below take into account the possibility that {B) ^ and 
(U) 7^ 0. The basic state therefore can be the outcome of a saturated small-scale 
dynamo, or it can result from the presence of a mean flux, as in the 2D case. It can 
of course be a mixture of both. 

Again, we are interested in the growth rate of long-wavelength perturbations, 
with wave vector k, which we consider to be of the form 



{b{x,t),u{x,t),p{x,t)) = {H{x,t),V{x,t),n{x,t)) ^k-x+p(k)t^ ^Al) 



where V, H and H vary on the same spatial and temporal scales as the basic state, 
and p{k) is the fe-dependent growth rate of the perturbations. The wavelength of 
the perturbations, 27r/|fc|, is assumed to be long in comparison with the spatial 
periodicity of the basic state. We now decompose the perturbations into mean and 
fluctuating parts, with the average defined by 



where ( ) denotes an average over the spatial and temporal periodicities of the basic 
state; {H), etc. are thus constants. Our aim is to employ a mean field approach 
in order to determine the growth rate p{k) of modes with non-zero mean magnetic 
fields, i.e. modes for which {H) ^ 0. 
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The evolution equations for the perturbations (cf. (2.7) - (2.9)) are 

(p + Rm-^k'^) H+{dt- Rm-^V"^) H = {V + ik) x (V x B + U x H) 

+ 2iRm-^k-\7H, (A3) 

{p + Re-^k^) V +{dt- Re-^V^) V = {V + ik) ■ {HB + BH - VU UV) 

- Vn - ikU + 2iRe-^k ■ VV, (A 4) 

V-H + ik-H = 0, V -V + ik-V = 0. (A 5) 

As in the 2D analysis, we use k = \k\ as a small parameter and expand H, V, II 
and the growth rate p in powers of k, so that, for example, 

H = Ho + Hi + ...Hn + ..., (A6) 

where i?„ is of n*'' order in the components of fc. Each Hn can be decomposed 
into its mean and fluctuating parts, Hn = {Hn) + hn, similarly, Vn = {Vn) +Vn- 
At zeroth order, equations (A 3) - (A 5) become 

PoHo + {dt - Rm^^V^) Hq = W x (Vq x B + U x Ho) , (A 7) 
PoVo + {dt - Re-'V^) Vo = -VHo + V • {HoB + BHq - VqU - UVq) , (A 8) 

V-Ho = 0, V-To = 0. (A9) 

On averaging these equations, we obtain 

Po{Ho)=Po{Vo} = 0. (A 10) 

For non-zero mean field solutions we therefore need to take po = 0. The correspond- 
ing equations for the fluctuations (cf. (2.15) - (2.17)) then read 

{dt - Rm-^V'^)ho = Vx{voxB + Uxho) + {Ho) ■ VU - {Vo) ■ VB, (A 11) 
{dt - Re-'^V'^)vo = -VHo + V • {hoB + Bho - vqU - Uvq) 

+ {Ho) ■VB-{Vo)-VU, (A 12) 

V-/io=0, V'-uo^O. (A 13) 

To first order, the equations for the mean variables (cf. (2.19) - (2.20)) give 

Pi{Ho) =ikx {vqX B + U X ho), (A 14) 

(Vo) = -ikUo + ik- {hoB + Bho - voU - Uvo) , (A 15) 

k-{Ho)=0, k-{Vo)=0. (A 16) 

It can be seen from equations (A 11) (A 12) that ho and Vq are subject to a linear 
forcing by both {Hq) and (Vq). Therefore, expressions (A 14) and (A 15) can be 
rewritten as 

Pi{Ho) =ikx (a^ • {Ho) + a" ■ {Vo)) , (A 17) 

Pi{Vo) = -ikilo + ik ■ (r^ • {Ho) + • (Vo)) , (A 18) 
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B 


Relevant matrices 





A = i 


-1.7667 -1.7684 
1.7684 1.7667 




0.033 0.386 

-0.386 -0.033 




0.1 


M = i 


" -1.311 -1.297 -0.376 0.405 " 
1.310 1.294 0.486 -0.295 
0.203 0.236 0.067 0.296 
0.273 0.229 -0.305 -0.038 





Table 4. Components of the matrices A., Q and Ai. for the AKA forcing and different 

values of B; see section 3. 



Rm 



16 



B 



0.001 



M 

0.988 
-0.988 











0.379 








i 


-0.019 






















0.815 










-0.202 






0.98 
-0.988 



0.25 



0.690 
-0.565 



1.0 





0.038 






1.228 














0.375 


-0.509 





0.495 












0.053 
-0.145 



128 



0.1 





-0.237 







1.140 



1.429 




1.233 
-1.0 



Table 5. Components of the matrices M. and A{B) for the MW+ forcing and different 

values ofB; see section 4- 



where the tensors a , or , T and F are constant and depend on the basic state 
U, B, and hence on the forcing F and on the Reynolds numbers. 

These four tensors need to be evaluated for an accurate determination of the 
growth rate of the perturbations to first order in fc. In theory, they can be calculated 
by solving equations (A 14) - (A 16). In practice however, these equations might 
have exponentially growing solutions independently of the forcing due to the mean 
magnetic and velocity fields (though if they are considered as long time solutions 
of an initial value problem this does not seem likely); we then may encounter the 
same problem as undermines the determination of the kinematic a-effect in the case 
when a small-scale dynamo is present (Cattaneo & Hughes 2009). 



Appendix B. Components of the matrices 

In tables 4 and 5 we detail the components of the matrices whose eigenvalues are 
given in the main text, for the AKA forcing and the MW-|- forcing respectively. 
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